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In this paper we investigate the chameleon effect in the different conformal frames of the Brans- 
Dicke theory. Given that, in the standard literature on the subject, the chameleon is described 
in the Einstein frame almost exclusively, here we pay special attention to the description of this 
effect in the Jordan and in the string frames. It is shown that, in general, terrestrial and solar 
system bounds on the mass of the BD scalar field, and bounds of cosmological origin, are difficult 
to reconcile at once through a single chameleon potential. We point out that, in a cosmological 
context, provided that the effective chameleon potential has a minimum within a region of constant 
density of matter, the Brans-Dicke theory transmutes into general relativity with a cosmological 
constant, in that region. This result, however, can be only locally valid. In cosmological settings 
de Sitter-general relativity is a global attractor of the Brans-Dicke theory only for the quadratic 
potential V(<j >) = M 2 <j> 2 , or for potentials that asymptote to M 2 cf> 2 . 

PACS numbers: 02.30.Jr, 02.30.Mv, 04.20.Fy, 04.60.Cf, 11.25.Wx 


I. INTRODUCTION 

The Brans-Dicke (BD) theory of gravity [lj represents 
the simplest modification of general relativity (GR): in 
addition to the 10 degrees which are associated with 
the metric tensor g^ v , a new scalar degree of freedom 
is also responsible for propagating gravity. This theory 
has been cornerstone for a better understanding of sev¬ 
eral other modifications of general relativity, such as the 
f(R) theories of gravity Q. In contrast to Einstein’s GR, 
the BD theory is not a fully geometrical theory of grav¬ 
ity, since, while one of the carriers of the gravitational 
field: the metric tensor, defines the metric properties 
of the spacetime, the scalar field <f>, which modifies the 
local strength of the gravitational interactions through 
the effective gravitational coupling G e s oc </> , is a non¬ 

geometric field. 

Although many aspects of BD theory have been well- 
explored in the past a a, other aspects have been 
cleared up just recently. Thanks to the chameleon ef¬ 
fect 0-113) f° r distance, it was just recently understood 
that the experimental bounds on the BD coupling param¬ 
eter w BD , which were set up through experiments in the 
solar system, might not apply in the large cosmological 
scales if consider BD theory with a potential. According 
to the chameleon effect, the effective mass of the scalar 
field 77i0, depends on the background energy density of 
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the environment: In the large cosmological scales where 
the background energy density is of the order of the crit- 
ical density p C rit ~ 10 -31 g/cm 3 , the effective mass is 
very small ~ H 0 ~ 10 -33 eV, so that the scalar field 
has impact in the cosmological dynamics. Meanwhile, in 
the solar system, where the averaged energy density of 
the environment is huge compared with p C rit, the effec¬ 
tive mass is large > 1 mm -1 {m<f, > 10 -3 eV), so 
that the Yukawa-like contribution of the scalar field to 
the gravitational interaction oc e~ m * r /r, is short-ranged, 
leading to an effective screening of the scalar field in the 
solar system. 

There is one aspect of the chameleon effect that has 
not yet been discussed in detail. If take a look at the ex¬ 
isting bibliography on this subject one immediately finds 
that this effect is almost exclusively described in the Ein¬ 
stein frame (EF), where the scalar field (the chameleon) 
is minimally coupled to the curvature scalar, at the cost 
of being non-minimally coupled to the matter sector of 
the action. A detailed description of this effect either in 
the Jordan frame (JF) or in the string frame (SF) is yet 
lacking. Since the chameleon effect is apparent in the 
density dependence of the dilaton’s mass, we think that 
the absence of appropriate discussion in the JF/SF, is 
due to the unconventional way which the self-interaction 
potential of the dilaton arises in the JF/SF Klein-Gordon 
(KG) equation, that governs its dinamics. 

In the present paper we shall try to fill the gap and the 
focus will be in the description of the chameleon effect in 
the Jordan and in the string frames of the BD theory with 
a potential (see sections ITO Mil HV1 and I VIII) . Although 
the very warped subject of the conformal transformations 
controversy will not be discussed here, the consequences 
for the chameleon effect of the most widespread view- 
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points on the (in)equivalence of the JF/SF and the EF, 
will be briefly discussed (see sections IXI1 IXII1 and IXIIID . 

In section HXl we will show that, in a cosmological con¬ 
text, provided that the effective chameleon potential has 
a minimum within a region of constant matter density, 
the GR-de Sitter solution can be, at most, either a lo¬ 
cal attractor or a saddle point of the Brans-Dicke theory 
within that region. In contrast, as it has been shown 
in OMzl by means of the tools of the dynamical sys¬ 
tems theory, in a cosmological setting the GR-de Sitter 
solution can be a global attractor of the BD theory ex¬ 
clusively for the quadratic potential: V{<j>) = M 2 * ^ 2 , or 
for any BD potential that asymptotes to the quadratic 
one, !/(</>)—>• M 2 (j) 2 (see the section IVIIII) . 1 

In section [VI] of this paper (see also the section IXl). we 
will show that the expectation that the chameleon effect 
might have had for the Brans- Dicke theory, specially for 
the relaxation of the stringent lower bounds on the value 
of the BD coupling constant w BD > 40000, is unjustified 
in general. It will be shown in sections [V] and IVIII (see 
the discussion in IXIVI) . that several BD models of im¬ 
portance for cosmology: those whose dynamics is driven 
by the quadratic potential oc </>, or by potentials 

that asymptote to the quadratic one, 2 do not develop the 
chameleon effect. This latter effect is primordial for the 
relaxation of the bounds on w BD , since, what one expects 
is that, thanks to the chameleon effect, the cosmologi¬ 
cal BD field is screened from terrestrial and solar system 
experimentation, even for w B d ~ 1. 


87tG = 1 /M| L = c = 1, are adopted. The field equations 
which are derived from m are the following: 


= 


W B D 


dii<t>d v (j) - ^g^ (< 9<^) 2 


9nv- 


V 


+i(V ll d v <l>-g liV V 2 <f>)+^T^\ ( 2 ) 


V ^=inkr»(^ v - 2V+1 2 TM )' 


( 3 ) 


where G^ = R— g^R/ 2, is the the Einstein’s tensor, 
V 2 = is the D’Alembertian operator, and 


rp(m) 

^ flU 


2 ^ (vTI^m) 

VW\ 99^ 


is the conserved stress-energy tensor of the matter de¬ 
grees of freedom 


WT™ = o. 

It is also convenient to rescale the BD scalar field and, 
consequently, the self-interaction potential: 


<t> = e v , V(4>) = e* U(<p), (4) 

so that, the action m is transformed into the string frame 
BD action: 


II. BASIC SETUP 

Here we assume the Brans- Dicke theory [l] with the 
potential, to dictate the dynamics of gravity and matter. 
In the Jordan frame it is depicted by the following action: 


S& = j d 4 xy/\i\^ {7? - w BD (a^) 2 - 2 u + 2 e~ v C m } .(5) 
The following motion equations are obtained from ([5]) : 


Sf F = y'd 4 a-v^l|^-^W) 2 -2U + 2£ m |, (1) 

where ( dq i) 2 = g^d^dv(j), V = V{<j>) is the scalar field 
self-interaction potential, w BD is the BD coupling param¬ 
eter, and C m = C m (x, d\, gpv) is the Lagrangian density 
of the matter degrees of freedom, collectively denoted 
by X- Unless the contrary is specified, the natural units 


G^v — (a 


1) 


d^dM- -g^(dip) 2 


\^) 2 + u(ip) 


-9^ 

2 


/iu f 


S7 2 y+{dvY = 


d v U -U+ 


3 + 


( 6 ) 

( 7 ) 


1 There are found several works on the de Sitter (inflationary) 
solutions within the scalar—tensor theory in the bibliography, in 
particular within the Brans-dicke theory, but just for illustration 
here we mention those in the references [181119t| . 

2 As shown in section rvTTTl these are the only BD cosmological 

models that have the ACDM model as a global attractor. 


where V 2 = , G^ = R^ - g^R/ 2, and, as 

before, is the (conserved) stress-energy tensor of 

the matter degrees of freedom: = 0. 

It is clear that the equations (HJ, J2J, ©, in the one 
hand, and ©, <[6|), dTJ) , in the other, transform into each 
other under the redefinitions in Eq. ©■ 
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III. THE KLEIN-GORDON EQUATION AND 
THE MASS OF THE SCALAR FIELD 

The mass (squared) of the BD scalar field can be com¬ 
puted with the help of the following equation [2(ij |: 

m\ = 3 + 2ta BD ~ ' ( 8 ) 

This mass is the one which is associated with a Yukawa- 
like term cf>(r) oc exp(— m^,r)/r, when the Klein-Gordon 
equation ([3]) is considered in the weak-field, slow-motion 
regime, and in the spherically symmetric case. The way 
in which this result is obtained can be found in the refer¬ 
ence [2(| in all details. However, for completeness of our 
exposition, here we shall explain the main reasoning line 
behind this result. 

In general for a scalar field which satisfies the standard 
KG equation 


V 2 </» = 5 0 U eff + 5, (9) 

where V e g = V e g{<j)) is the effective self-interaction po¬ 
tential of the scalar field <j>, and S is a source term (S 
does not depend on <f>), the effective mass squared of the 
scalar field is defined by m| = The problem with 

this definition is that the BD scalar field does not satisfy 
the usual KG equation, but the one in Eq. ©, where 
the self-interaction potential of the BD field appears in 
an unconventional way. In order to fix this problem, one 
notices that, if introduce the effective potential (2C3J 


V e g{(t>) = 


3 T 2w B 


4>V{^>) - 3 / d(j)V{<j)) 


( 10 ) 


so that 

d^V eS = 2 [^V{4>) - 2!/(</>)], 

O “T 

then, the KG equation m can be rewritten in the more 
conventional way: 


V 2 0 = d^Ves + »,i (11) 

O T" 

where the second term in the right-hand side (RHS) of 
this equation, is the source term which does not depend 
explicitly on the 0—field. 


A. Perturbations around the minimum of the 
effective potential 

Usually the details of the derivation that leads to the 
appearance of an effective mass of the scalar field are 


not given. Here we want to include these details, since 
we want to make clear that, what one usually calls as an 
effective mass, is a concept that is linked with the oscilla¬ 
tions of the field around the minimum of the effective po¬ 
tential, which propagate in spacetime. These oscillations, 
or excitations, are the ones that carry energy-momentum 
and, if required, may be quantized. 

For simplicity consider the vacuum case = 0 of 

Eq. (ED. Let us assume next that V e g is a minimum at 
some 4>* ■ Given that the equation (ED is non-linear, one 
may consider small deviations around </>*: </> = <^>* + S(f> 
{5(f) -C 1), then, up to terms linear in the deviation, one 
gets: 

fyVefffa) ~ d<t>Ve+ dlV e s(4>*) $4> + ■■■ = TO*<I</>, 

where to* = d^Ves (<(>*), is the effective mass of the scalar 
field perturbations. 

Working in a flat Minkowski background (in spherical 
coordinates) 

ds 2 = —dt 2 + dr 2 + r 2 (dO 2 + sin 2 ddcj) 2 ) , 

which amounts to ignoring the curvature effects and the 
backreaction of the scalar field perturbations on the met¬ 
ric, and imposing separation of variables, the perturba¬ 
tions 


6<p = 5ip(t,r) = J2 e “”Vn(r), 

n 

where w n is the angular frequency of the oscillations of 
the n-th excitation of the field, obey the Helmholtz equa¬ 
tion: 


d 2 i> n 

dr 2 


, 2 d,ip n 2 

d--k krj'i/'n — 0, 

r ctr 


( 12 ) 


where = w 2 — to*, is the wave-number squared. Eq. 
El is solved by the spherical waves: 


pi| k nb 

ipn{r) = c n - 


where the C n - s are integration constants. For |w„| < to*, 
the Eq. (fl2l) has the Yukawa-type solution: 


— — uj 2 r 

ipn(r) = C n -. (13) 

r 

It is understood that the modes with the lowest en¬ 
ergies E n = hu} n = u} n <C to*, are the ones which are 
more easily excited in the small oscillations approxima¬ 
tion, and, hence, are the prevailing ones. In particular, 
wo -C to*, so that 
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-m„ r 

ip(r) ~ ip 0 (r) = C —-—. (14) 

These lowest order modes are the ones with the shortest 
effective Compton wave length A* ~ to" 1 , and are the 
ones which decide the range of the Yukawa-type interac¬ 
tion, i. e., these are the ones which decide the effective 
screening of the </>-field. 

In what follows, as it is usually done in the related bib¬ 
liography, when we talk about perturbations of the scalar 
field around the minimum of the effective potential, we 
ignore the oscillations in time by assuming the static situ¬ 
ation, etc. This amounts, effectively, to ignoring all of the 
higher-energy excitations of the field. This assumption 
bears no consequences for the qualitative analysis. How¬ 
ever, once a friction term oc <f> arises (the dot accounts 
for the time derivative), for instance, in a cosmological 
context, the oscillations of the field in time, around the 
minimum, are necessarily to be considered (see section 

m- 


B. The field—theoretical mass 

The above analysis suggests that one may introduce a 
general definition of the mass (squared) of the BD scalar 
field ml^ml= <9^T4t (</>*), 

ml = dlV eS {4>) = x— \ - y>dlV{<t>) - d 0 V(c£)] , 

O “h 

which is just Eq. © at the beginning of this section. 

What if the effective potential V e g has no minimums 
at all? The quartic potential = A</> 4 , for instance, 

has a minimum at <j> = 0, however, the corresponding 
effective potential: 


Kff(^) = 


4A </> 5 

5(3 + 2w BD ) ’ 


(15) 


has no minimums. In this case, our understanding of 
what an effective mass means, might have no meaning 
at all. In particular, the screened Coulomb-type poten¬ 
tial (the mentioned Yukawa-like solution), being the most 
relevant physical manifestation of a massive propagator, 
might not arise. The corresponding “mass” in Eq. © 
would be just a useful field theoretical construction with 
the dimensions of mass, no more. 

In spite of this, following the most widespread point of 
view on the correct definition of the mass of the BD field 
in the Jordan frame (2Cj, here we shall consider that, in 
general, i. e., even away from the minimum of the ef¬ 
fective potential, the field parameter m| given by Eq. 


©, represents the effective mass (squared) of the scalar 
field. In this paper, in order to differentiate the men¬ 
tioned field theoretic parameter from an actual effective 
mass (the one with consequences for fifth force experi¬ 
ments) , we shall call the latter as “effective mass”, while 
the former as “effective field-theoretical mass” (see the 
related discussion in section IXIII CD . 


IV. THE CHAMELEON MASS 

Due to the chameleon mechanism, the screening effect 
may arise even if the effective potential V e g does not de¬ 
velop minimums. As we shall see, all what one needs is to 
include the source term in the RHS of the KG equation 
© , within a redefined effective potential, which we shall 
call effective chameleonic potential: 


V c h = Keff + (j>S =► V 2 </> = d^V ch = d^V eS + S. (16) 

Usually the chameleon effect is discussed, exclusively, 
in the Einstein frame formulation of the Brans-Dicke 
theory, where the scalar field couples directly with the 
matter degrees of freedom dUH. Due to the non—trivial 
way which the self interaction potential enters in the KG 
equation, the discussion of the chameleon effect in the 
Jordan frame formulation of BD theory seems more ob¬ 
scure than in the Einstein frame. 

Here we shall show that, regardless of the unconven¬ 
tional form of the potential in the BD-KG equation © , 
the chameleon effect can be discussed in the Jordan frame 
as well, if in Eq. © introduce the following definition of 
the effective chameleon potential: 


U ch (^) = V eS (cf>) 


1 


3 J- 2w b 


■<^>T (m) = 


3 + 2 id# 




,( 17 ) 


so that 


= - 2V + i!■'")) , (18) 

coincides with the RHS of the KG Eq. ©, and the lat¬ 
ter can be written in the form of the conventional KG 
equation without a source: V 2 </> = d^V c h. 

As in the standard case, the effective mass squared: 
m l, = <9|l4h (</>*), may be defined for the small pertur¬ 
bations of the BD scalar field around the minimum 
of the chameleon potential I4h- Actually, under the as¬ 
sumption of spherical symmetry, given that V c h is a min¬ 
imum at some </>*, if follow the procedure explained in 
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the former section Hill in the weak-field and low-velocity 
limit (basically the case when the curvature effects and 
the backreaction on the metric are ignored): 


where 


d 2 ( 64 >) 
dr 2 


2 d(6<t>) 

r dr 




m l, = d l V ch {(/>*), 

is the effective mass of the perturbations around the min¬ 
imum of the chameleon potential. 

We recall that, as shown in the section Ull Ai although 
when dealing with the chameleon effect we care only 
about the spatial deviation about the minimum <50(r), 
as a means to linearize the Klein -Gordon equation, in 
general these deviations are also time-dependent, so that 
what we have are time-dependent perturbations around 
the minimum of the chameleon potential. These may be 
viewed as periodic oscillations of the BD field about the 
minimum, and the resulting effective mass can be inter¬ 
preted as the mass of the corresponding scalar excitations 
propagating in a flat background. 

Solving for the above Helmholtz equation, one has for 
<j>(r) = 0* + 8<j>(r), the following solution: 

P ~ m 4,, r p m <K»~ 

</>(r) = </>*+ Ci -b C 2 -, 

r r 

where C\ and C 2 are integration constants which we can 
determine through the boundary conditions. If assume, 
for instance, that 0(r) tends asymptotically to a constant 
value (j) 00 ■ 

lim 0(r) = 4 >oo => <t>(r) = (j)00 + C\ -. 

r—too J‘ 

Depending on the physical situation at hand, other 
boundary conditions are required in order to fix the re¬ 
maining constants C\ and Gi- For instance, if assume 
regularity of the solution at the origin r = 0, then 
Ci = -C 2 , 

At \ _ a , n smh(m^r) 

(pyr) — 0 * + C'2-, 

r 

so that 0(0) = 0* + C 2 TO*, etc. 

As we shall see, the interesting thing here is that the 
effective chameleon mass m ^ (p), is a function of 

the surrounding density p. This property of the effective 
mass of the BD scalar field perturbations is what is called, 
primarily, as the chameleon effect. 


A. On the matter density 


argument of the effective mass: m^„(p), is the density 
measured by co-moving observers (with four-velocity <5g) 

in the JF-BD theory: p = T^SqSq. This is, besides, 
the density of matter that is conserved in the JF (also in 
the SF) formulation of the Brans Dicke theory. 

This is to be contrasted with the original chameleon 
effect of Ref. 0, where the physically meaningful mat¬ 
ter density pi is not the one measured by co-moving 
observers in the EF (this is not the conserved one in 
this frame), neither in the conformal one, but a den¬ 
sity which does not depend on the dilaton. Actually, 
in 0 the matter fields couple to the conformal metric 
gfy = exp(2/3j0/M PL )g M1/ , while the density of the non- 
relativistic fluid measured by EF co-moving observers is 
denoted by pi. It is assumed that what matters is the 0- 
independent density p,; = pi exp(3/3j0/M PL ), which is the 
one conserved in the EF. While this choice may not be 
unique, in the Jordan frame of the Brans-Dicke theory 
(the same for the SF), one does not have this ambigu¬ 
ity: the matter density measured by JF(SF) co-moving 
observers p, is the one conserved in the Jordan/string 
frames and, additionally, does not depend on the BD- 
field. 


V. THE BD CHAMELEON: EXAMPLES 

In this section we shall to illustrate how the chameleon 
effect arises in the Jordan frame of the Brans- Dicke the¬ 
ory, by exploring a pair of examples. 


A. The quartic potential as an example 


In order to illustrate our reasoning 
example with the quartic potential [6| 


let us choose the 


F(0) = A0 4 , (19) 

where we assume that the free parameter A > 0, is a 
non-negative constant. In this case, as said, the effective 
potential (fl5l) : I4ff(0) oc 0, does not develop minimums. 
Yet the corresponding chameleon potential 03: 


F ch (0) = 


4 

3 + 2w B d 




4 

3 + 2w B d 




( 20 ) 


An interesting aspect of the BD chameleon effect is 
related with the fact that the density of matter p in the 


where we have assumed a homogeneous, pressureless dust 
background: = — p, can have a minimum if w BD > 
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—3/2. Actually, at the value 


*• = («) 


1/4 


the derivatives of the above chameleon potential 

16A /on 3 / 4 


d</>Vch(<f>*) = 0 , (<£*) = 


3 A 2c0b 


(JLY 

V4A/ ' 


Hence, provided that w BD > —3/2, since <9^V/h((/*) > 0, 
the chameleon potential V c h in Eq. m , is a minimum at 
</>„. In this case we can identify a physically meaningful 
effective mass of the BD field: 


m l, = d|Kh(0*) 


16A /_P_\ 3/4 

3 v 4A / 


( 21 ) 


We want to underline that, in general, p can be a func¬ 
tion of the spacetime point p = p(x), however, in most 
applications the function p(x) is assumed piece-wise con¬ 
stant. For instance, one may imagine an spherical spatial 
region of radius R , filled with a static fluid with homoge¬ 
neous and isotropic constant density po, and surrounded 
by a fluid with a different (also homogeneous) constant 
density poo , so that: 


p[r) 


Po if r < A; 
Poo if r » R. 


In such a case the effective mass of the BD scalar 
field would be one for modes propagating inside the 
spherical region and another different value m^, 
for scalar modes propagating outside of (far from) the 
spherical region: 

W0 O = yJd 2 V ch ((j ) 0 ), = yJd 2 V c hO/oo). 

For the quartic potential, in particular, one would have 
that: 


_ 2(4A) 1 / 8 pl /S _ 2(4A) 1 / 8 p^j S 

77100 V3 + 2w B d ’ m0 “ V3 + 2w bd ’ 

inside and outside of the spherical region of radius R, 
respectively. 

In case the the gravitational configuration of matter 
were given by a point-dependent density profile p = p(x), 
such as, for instance, in a cosmological context where 
p = p(t) ( t is the cosmic time), the effective chameleon 
mass were point-dependent as well. However, as 
it is well known, the masses of point particles in the 
JF/SF formulations of the BD theory, are constants by 
definition. Otherwise, these particles would not follow 
geodesics of the JF/SF metric. In general, coexistence of 


particles of constant mass and particle excitations with 
point-dependent mass, bring about problems with the 
equivalence principle. Besides, if p = p{x), then, the re¬ 
sulting effective mass will be a field-theoretical construc¬ 
tion which, as it will be shown in section IXIII Cl has an 
anomalous behavior under the conformal transformation 
of the metric. 

In order to evade any possible discussion on the equiv¬ 
alence principle, or on the anomalous behavior of the 
effective field-theoretical chameleon mass under the con¬ 
formal transformations of the metric, here we adopt the 
most widespread handling of the chameleon effect, and 
we assume that the density of matter has piece-wise con¬ 
stant profile in the sense explained above. In section m 
we will discuss more on this delicate subject. 


B. The quadratic monomial: the massless BD held 


One peculiar note about the effective chameleon poten¬ 
tials: if look at equation (El), one sees that the quadratic 
potential V{(j)) = M 2 cj) 2 , plays a singular role. Actually, 
if substitute this potential into ED. one obtains that the 
resulting chameleon potential 


Vch(/>) 


<; bT ( m ) 

3 + 2w bd 


does not have a minimum. Besides, the second derivative 
vanishes: d^Vch = 0. This means that the quadratic 
potential does not generate the chameleon effect. The 
same is true for any potential that asymptotes to 4> , for 
instance V((f >) oc cosh(A0) — 1. Even the effective field- 
theoretical mass squared ©, vanishes for the quadratic 
monomial. 

As a consequence, assuming the ’mexican hat’ poten¬ 
tial 


V{4>) = —M 2 (j) 2 + A 0 4 , 

the corresponding effective and chameleon potentials 
m, and ED, will be given by: 

4A </> 5 _ 4 

eff 5(3 + 2 w bd ) ’ ch 3 + 2 w bd 

respectively, which coincide with the effective and 
chameleon potentials for the quartic monomial (1TD1) . This 
means that the quadratic monomial does not contribute 
neither to the effective nor to the chameleon potentials. 

In the case of a standard scalar field a whose dy¬ 
namics is governed by the usual Kelin-Gordon equation 
V 2 er = d a V(a), the quadractic monomial V(a) oc cr 2 , is 
also a singular potential in the sense that, it is the only 
potential for which the KG equation is a linear differential 
equation, i. e., the superposition principle is satisfied. 
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VI. THE QUARTIC POTENTIAL: ESTIMATES 


Notice that, similar to the chameleon effect arising in 
the Einstein frame of the BD theory 0, the mass squared 
of the BD field given by Eq. e., the Jordan frame 

mass - the one that determines the range of the Yukawa- 
like correction 1 211 1 - depends on the background energy 
density oc p 3 ' 8 . As it can be seen, this dependence 
of the mass of the scalar field on the ambient energy 
density improves the one in Ref. [G]: oc p 1 / 3 , just by 

a fraction. 

In order to make estimates, let us write 


4 5/8 A l/8 

TO0 * = V3 + 2w bd 



3/8 




or in 


“user-friendly” units (using the terminology of 


0 ): 


10A 1/8 / , \3/8 

i NitH. t*/" 1 ) ■ (22) 

Let us assume that the scalar field is immersed in the 
earth atmosphere with mean density p atm ps 10 -3 g/cm 3 , 
then, provided that the millimeter range screening 0: 
( m at in )-i ^ 4 mm, is undertaken, from Eq. (l22l) it fol¬ 
lows that 


w BD w 1.6 A 1/4 - 1.5, 

so that, if consider, for instance, that A ~ 1, one gets 
that w B d ~ 0.1 can be of order unity or smaller. What 
this means is that the BD theory may describe the grav¬ 
itational phenomena with a coupling constant of order 
unity and, yet, the chameleon potential EU1) may effec¬ 
tively screen the BD field from experiments that look for 
violation of the Newton’s law, for distances above the 
millimeter. 

The next question is whether the above potential can 
be a good candidate for cosmology as well. The ratio of 
the mass of the BD field measured at large cosmologi¬ 
cal scales, to the scalar field mass estimated in earth’s 
atmosphere: 


rrij 


3/8 


3 x 10 


-n 


(23) 


where we have taken into account that the critical energy 
density of the universe p crlt ~ 10~ 31 g/cm 3 . If consider 
the millimeter-range screening above, ss 1mm 

=>• toS 3 10~ 4 eV, then the estimated mass of the 
cosmological BD scalar field 

rn£ sm «3x 10 —11 TO^ 4m » 3 x 10” 15 eV, 


is by some 18 orders of magnitude heavier than the ex¬ 
pected value TOS° sm ~ H 0 ~ 10 -33 eV. Hence, if as¬ 
sume that the BD scalar field with a fixed potential 
V(4>) = A/> 4 , is effectively screened from solar system ex¬ 
perimentation, the BD field would not have cosmological 
implications. 

The “reconciliation” between terrestrial and cosmolog¬ 
ical bounds, at once, can be achieved by power law po¬ 
tentials leading to BD chameleon mass, oc ( p ) k / 2 , 
with the power k ss 29/14 ss 2.071, or higher. Of course, 
the reconciliation is natural if, for instance, oc exp p. 


A. Thin—shell effect 


Our estimates above are unsatisfactory in many as¬ 
pects. First of all, a lot of simplification has been made 
for sake of transparency of our analysis. For instance, 
the well known thin-shell effect 0 , which arises due to 
the non-linearity of the BD scalar field, and which is sig¬ 
nificant for large bodies, has not been considered in our 
analysis. 3 Nevertheless, even if take into account the 
thin-shell effect, the physical implications of the huge 
difference between the cosmic and terrestrial mass scales: 
TO cos m / TO atm ^ 10 —11 , can not be erased by the (thin 
shell mediated) weakening of the effective coupling of 
the BD field to the surrounding matter. Actually, the 
additional contribution of the chameleon BD field to the 
Newtonian gravitational potential energy of a given mass 
Mb, is expressed by 


p~ i / x'ett 

A U* N oc —-—, 

where /3* s is the effective coupling of the chameleon field 
to the surrounding matter, and A e ff = m^ 1 , is its effective 
Compton length. We have that 


A atm _ 10 18 A atm 5 (24) 

where A^* 111 ~ 1mm, is the Compton length of the 
chameleon field which is consistent with the experiments 
on fifth-force, while A|/ m is the effective range of the 
scalar field mediated interaction, computed with the po¬ 
tential V(<j>) oc q f> 4 , under the assumption that the cos¬ 
mological bound A^/ sm ~ 10 26 m, is met: 

A^‘ m ~ 1CT n A£ sm ~ 10 15 m = 10 18 mm. 


3 For a detailed exposition of the thin-shell effect we recommend 

Q (See also the appendix IX VII) . 













Then, requiring that 


A U* N (PtA 2 e^r 
AUn \/3efiJ e r / x tT 

for the given potential V(4>) oc 0 4 , the expected weak¬ 
ening of the effective coupling of the chameleon BD field 
to the surrounding matter, is an unnaturally large effect: 
(/3* ff ) 2 ~ exp (—10 20 ) /3 2 ff , where we have assumed that 
/3 e ff ~ 1. In order to obtain the above estimate, we have 
arbitrarily set the distance from the source of gravity 
r ss 10 2 A^ m , and the Eq. (HH1) has been considered. 

The above results are true, in general, for power-law 
potentials of arbitrary power: V(4>) oc <fr. In this lat¬ 
ter case, for the mass squared of the BD field, one gets: 
to 2 oc /o^ A—1 )/ A , where, as A —> oo, (A — 1)/A —> 1. 
This means that the latter power can never exceed unity 
(A— 1)/A< 1. Recall that, a necessary requirement, when 
the power- law potential V (</>) oc <fr is allowed to explain 
cosmological and terrestrial bounds at once, amounts to: 
(A - 1)/A > 2.071. 

Our conclusion is that, in general, terrestrial and solar 
system bounds on the mass of the BD scalar field, and 
bounds of cosmological origin, are difficult to reconcile 
through a single chameleon potential. 


VII. THE CHAMELEONIC MASS OF THE 
DILATON 

In several situations of interest, in particular when the 
focus is in the low-energy effective string theory, it is 
useful to turn to the field variables p, and U = U(p), 
in Eq. O- This choice singles out the string frame for¬ 
mulation of this theory. In terms of these variables the 
effective mass squared of the dilaton, i. e., the one which 
is obtained from Eq. © by substitution of ©, is written 
as: 


[ d l U &) ~ U{ip)\ . (25) 

However, as we will immediately show, the above expres¬ 
sion for the effective mass of the dilaton, has to be taken 
with caution. 

If follow, step by step, the procedure in sections mu 
and hyi then, the KG equation © can be rewritten as: 


VV + (dp) 2 = 5 v C/ ch , (26) 

where the dilatonic chameleon potential U c h = U c h(p), is 
defined as: 


U ch = 


3 T 2wb 


U- dpU - - e~ v T (m) 


(27) 


so that 

d v U - U + i , 

coincides with the RHS of Eq. 0. 

Assuming that U c h is a minimum at some p* , small de¬ 
viations about p*: p = </?*+ Sp (dp <C 1), d^Sp —> d^p 
=> (d5p) 2 ~ 0(5p 2 ), so that, the second term in the left- 
hand side (LHS) of Eq. (EH1) . may be ignored, and one 
ends up with the standard KG equation for the fluctua¬ 
tions of the dilaton: S7 2 8p = d v U c h. Since U c h develops a 
minimum at <p*, then, one may define the effective mass 
(squared) of the dilaton: m 2 ^ = d 2 U c h(p *), which is the 
one that determines the range of the Yukawa interaction 
in the string frame: 


dtnUcYl — 


3 -f- 2 Wg 


g —m v ,r 

8p(r) oc -. 

r 

Following the procedure exposed in the former sections 
eh and hyi a straightforward generalization 

to 2 (<P) = d 2 U ch (p) —t m 2 ^ = d 2 U c h(p*), 

leads to the following definition of the effective field- 
theoretical mass of the dilaton in the string frame: 


'■M = 


3 + 2c Uq 


dlU-d^U- 


■ (28) 


It is not Eq. (l25l) . but the latter equation ([28]) . the 
one to be contrasted with the effective field-theoretical 
mass parameter © of the BD scalar field in the Jordan 
frame. Nevertheless, since, as assumed, t/ c h develops a 
minimum at some p *, then 


dipUchip*) — 0 

d v U(p*) = U(p*)- l -e-*'T( m \ 

so that, 


d 2 U ch (p*) = 


3 + 2cj b 


d%,U(p*) - d v U(p*) - ^——T^ rn ' ) 


3 + 2cj b 


K%)-%)], (29) 


where the last equation above coincides with Eq. (l25l) . 
when given quantities are evaluated at the minimum 
of the chameleon potential (IT71) . 
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A. Chameleonic dilaton 

For sake of illustration, let us assume, as before, a pres- 
sureless, homogeneous and isotropic background with en¬ 
ergy density p. Consider an exponential potential 


U(ip)=M 2 e x,p , (30) 

where M and A are free constants. Notice that, under the 
replacements </> —> e v , M 4 ~ n M 2 , and n — 1 —> A, the 
latter potential in Eq. (1301) is mapped into the potential 

V(0) = M 4 -”(/) n , 

in terms of the JF-BD scalar field <f>. The particular 
value A = 3 (n = 4), corresponds to the quartic potential 
discussed above. According to Eq. (1271) . the resulting 
chameleon potential looks like: 


U cb (ip) = 


2M 2 (A - 1) 
A(3 + 2wbd) 



A pe~ v 
2M 2 (A — 1) 


(31) 


where, as before, we have assumed a homogeneous 
pressureless fluid with trace of the stress-energy tensor 
T( m ) = -p. 

For A > 1, the above potential is a minimum at 4 


ip* = In 


P 

2M 2 (A — 1) 


(32) 


where the effective mass squared of the fluctuations of 
the chameleonic dilaton is given by to 2 ^ = d 2 U c h(p*)'- 


(A + l) [2M 2 (A-1)] a+1 _a_ 
3 + 2w B d 


(33) 


For the particular value A = 3 which, as mentioned, 
corresponds to the quartic potential in the Jordan frame 
of BD theory: V((f>) oc </> 4 , the chameleonic mass of the 
dilaton goes like oc p 3 / 8 , which is consistent with 
the results of section IVll as one should expect. 

As explained in former sections, an interesting thing 
about the mass of the dilaton is related with the fact that, 
for the specific exponential potential U(ip) oc e v , which, 
in terms of the BD field <j>, corresponds to the quadratic 
monomial potential V(4>) oc <f) 2 , the chameleonic dilaton 
is massless. Actually, in this case, in Eq. (1301) one have 


4 For A < 1, the chameleon potential ED does not develop mini- 

mums at all. 


to set A = 1. When this value of the free parameter is 
substituted into Eq. (l33l) one gets m 2 ^ = 0. This con¬ 
clusion is correct even if consider the effective potential 
m- As shown in section IV Bl in the Jordan frame of 
the BD theory, where the effective mass squared of the 
BD field is given by Eq. ©, for the quadratic mono¬ 
mial V((f>) = M 2 (j) 2 / 2, which corresponds to the expo¬ 
nential potential for the dilaton U(ip) = M 2 exp(<p)/2, 
since d^V = M 2 4> , and d 2 V = M 2 , the mass squared © 
vanishes: 


= 3+ ^ bd [<t>9lV(<j>) - d^V((/))] = 0 . 

Hence, the potential 

U(ip) oc e v <(=> V((j>) oc (j) 2 , 

does not provides an effective mass for the BD scalar 
field. 

This means that, for the exponential potential m 
with the specific value of the free parameter A = 1 (or 
for the quadratic monomial if working with the Jordan 
frame BD variables), the chameleon effect does not work, 
and the dilaton can not be screened from solar systems 
experiments. This looks like a bad news for the dilaton 
since, as shown in Ref. 01 (see the next section), 
only for the exponential potential U(tp) = M 2 expip, or 
for potentials which asymptote to M 2 exp ip, the BD the¬ 
ory has the ACDM solution - also called as concordance 
model [H HI - as a global attractor. 


VIII. WHAT DO THE DYNAMICAL SYSTEMS 
HAVE TO SAY ABOUT THE BRANS DICKE 
CHAMELEON? 


In this section we shall consider Friedmann-Robertson- 
Walker (FRW) spacetimes with flat spatial sections for 
which the line-element takes the simple form: 

ds 2 = —dt 2 + a 2 {t)8ijdx l dx : ’, i,j = 1,2,3. 

We assume the matter content of the Universe in the form 
of a cosmological perfect fluid, which is characterized by 
the following state equation p m = w m p m , relating the 
barotropic pressure p m and the energy density p m of the 
fluid, where w m is the so called equation of state (EOS) 
parameter. Under these assumptions the cosmological 
equations which are derived from ([5)1, are written as it 
follows: 
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3 H 2 = ^ip 2 - 3 Hip + U + e" 


jV Wed . 2 


2 Hip 


d v U-U 
3 + 2w bd 


2 + Wed (1 + Wm) _ 


3 T 2 a; i, 


^ ^ Pm ; 


W 

Prr 


2, Hip + ip 2 = 2 

+ 3iJ (m 


U — dtnU 1 — 3 w„ 


3 + 2w bd 

+ 1) Pm = 0, 


3 T 2cl?e 


° h'm •> 


(34) 


where H = a/a is the Hubble parameter. 

Due to the complexity of the system of non-linear 
second-order differential equations (1M1) . it is a very dif¬ 
ficult (and perhaps unsuccessful) task to find exact so¬ 
lutions. Even when an analytic solution can be found 
it will not be unique but just one in a large set of them. 
This is in addition to the problem of the stability of given 
solutions. In this case the dynamical systems tools come 
to our rescue. The dynamical systems theory provides 
powerful tools which are commonly used in cosmology to 
extract essential information on the dynamical properties 
of a variety of cosmological models, in particular, those 
models where the scalar field plays a role [lfl I24-I43j . 

In general, when one deals with BD cosmological mod¬ 
els it is customary to choose the following variables 0: 


_ ip ip’ _ VU _ d v U 

X ~ V6H ~ vr y ~ y/ZH' e U 


(35) 


where the tilde means derivative with respect to the vari¬ 
able r = In a, i. e., to the number of e-foldings. As a 
matter of fact x and y in Eq. (1351) . are the same vari¬ 
ables which are usually considered in similar dynamical 
systems studies of FRW cosmology, within the frame of 
Einstein’s general relativity with a scalar field matter 
source [24[. In terms of the above variables the Fried¬ 
mann constraint can be written as 0 


ft™ ~ =1 + ^ x ~ Wbd x2 ~~ y2 ~ °' 

Notice that one might define a dimensionless potential 
energy density and an “effective kinetic” energy density 

ft u = = y 2 , ft e K = x (ubbX - Tfi) , (37) 

respectively, so that the Friedmann constraint can be re¬ 
written in the following compact form: 

+ fi/y + ft^f = 1. 


The definition for the dimensionless effective kinetic 
energy density has not the same meaning as in GR 
with a scalar field. It may be a negative quantity with¬ 
out challenging the known laws of physics. Besides, since 
there is not restriction on the sign of then, it might 
happen that Dj/ = U/3H 2 > 1. This is due to the fact 
that the dilaton field in the BD theory is not a standard 
matter field but it is a part of the gravitational field it¬ 
self. This effective (dimensionless) kinetic energy density 
vanishes whenever: 

V6 6 6 

x = - => ip =- H => ip =- in a, 

Wbd W B d W B d 

or if: £ = 0 => ip = 0 => ip = const., which, provided 
that the matter fluid is cold dark matter, corresponds to 
the GR-de Sitter universe, i. e., to the ACDM model. 5 


A. The dynamical system 

Our goal will be to write the resulting system of cos¬ 
mological equations (1341) . in the form of a system of 
autonomous ordinary differential equations (ODE-s) in 
terms of the variables x, y , £, of some phase space. We 
obtain the following autonomous system of ODE-s: 


x' = — 3x (l + VGx — w BD £ 2 ) + 


x \/2/3 2 . 


3 + 2w B d 

+ [2 + w BD (l + r^m)] x 


3y £ 


l-3w 


qo eff 

05 L m > 


y =y 


3x ^w BD r — 


3 -f- 2oj,. 

£ + 3 


3 y 2 tj 
3 T 2u> b 


76 

2 + W BD (1 + Wm) ooeff 

3 + 2w bd 


7 = -7mi-0 2 (f-i), 


(38) 


5 It is known that, but for some anomalies in the power spectrum of 
the cosmic microwave background [44l| , at the present stage of the 
cosmic evolution, any cosmological model has to approach to the 
so called concordance or ACDM model [22|| . The mathematical 
basis for the latter is the Einstein—Hilbert action plus a matter 
action piece: 

'Sacdm = ~2 J d^x\/\g\ (R — 2A) + J d 4 xyj\g\Ccvu-, 

where £cdm is the Lagrangian density of (pressureless) cold dark 
matter (CDM). The ACDM action above can be obtained from 
provided that the dilaton acquires some expectation value 
ipo, so that 

= —, U(ipo) = A, £cdm £m • 
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where is given by Eq. HMD, and it is assumed that 
T = Ud 2 U/(d^U ) 2 can be written as a function of £ [29| : 
r = T(£). Hence, the properties of the dynamical system 
(1381) are highly dependent on the specific functional form 
of the potential U = U(ip). 

B. Vacuum Brans—Dicke cosmology 

A significant simplification of the dynamical equa¬ 
tions is achieved when matter degrees of freedom are 
not considered. In this case, since = 0 => y 2 = 
1 + \/6x — w BD x 2 , then the system of ODE-s (1551) simpli¬ 
fies to a plane-autonomous system of ODE-s: 

= (~ 3x + 3 " 3 +^ *) i 1 + ~ UB ° X2 ) ’ 

^ = - V6ar (1 - O 2 (E - 1) - (39) 

In the present case one has 

= 7^2 = y 2 = 1 + 'f§x - u BD x 2 , 

nf = x (w BD x - y/6j => nf+Qu = 1, (40) 

where we recall that the definition of the effective (di¬ 
mensionless) kinetic energy density has not the same 
meaning as in GR with scalar field matter, and it may be, 
even, a negative quantity. Here we consider non-negative 
self-interaction potentials U(<p) > 0, so that the dimen¬ 
sionless potential energy density flu = y 2 , is restricted 
to be always non-negative: flu = 1 + \/6x — ui B D x 2 > 0. 
Otherwise, y 2 < 0, and the phase-plane would be a com¬ 
plex plane. Besides, we shall be interested in expanding 
cosmological solutions exclusively (H > 0), so that y > 0. 
Because of this the variable x is bounded to take values 
within the following interval: 

_ _ /3 (l ± ylT2 w bd /3^ 

a - — x — Q +> Q± = V 2 (- TTn - )' <41) 

This means that the phase space for the vacuum Brans- 
Dicke theory \l/ vac can be defined as follows: 

^vac = {(*, 0 : a- <x < a + } , (42) 

where the bounds on the variable £ - if any - are set by 
the concrete form of the self-interaction potential. 

There are found four dilatonic equilibrium points Pi : 
( Xi , £,;), in the phase space corresponding to the dynami¬ 
cal system m- One of them is the GR-de Sitter phase: 


PdS : (0,0) => x = 0 => tp = (p 0 , and 
y 2 = 1 => 3 H 2 = U = const., 

which, since H given by 

Jp = ( 2 V 6 - 3w bd x) X 

3(l + V&x — uj bd x 2 ) £ 

3 + 2w bd ’ ( ’ 

vanishes, then H = 0 => H = H$, corresponds to accel¬ 
erated expansion q = — 1 — H/H 2 = — 1. Besides, 

f = 0 ^ = 1 =► U(ip) oc e v , 

i. e., this point exists for this specific exponential po¬ 
tential exclusively. The eigenvalues of the linearization 
matrix around (0, 0) are: Ai = —3, A 2 = 0. This means 
that (0, 0) is a non-hyperbolic point. 

We found, also, another de Sitter solution: q = — 1 => 
H = 0, which is associated with scaling of the effective 
kinetic and potential energies of the dilaton: 

pBD . f 1 . Of = 6 + 5 o; bd 

dS \-\/6(l + Wed) ' / tlu 12 + 17w B D + 6w| D ’ 

Ai = —(4 + 3 w bd )/(1 + w BD ), A 2 = 0, where, as before, 
Ai and A 2 are the eigenvalues of the linearization matrix 
around the critical point. We call this as BD-de Sitter 
critical point to differentiate it from the GR-de Sitter 
point. 

In order to make clear what the difference is between 
both de Sitter solutions, 6 let us note that the Fried¬ 
mann constraint (l36l) . evaluated at the BD-de Sitter 
point above, can be written as 

e^p m = 3H 2 + - 

i. e., e~ v p m = const. This means that the weaken¬ 
ing/strengthening of the effective gravitational coupling 
(G e fi oc e - ^) is accompanied by a compensating grow¬ 
ing/decreasing property of the energy density of matter 
p m oc e v , which leads to an exponential rate of expan¬ 
sion a(t) oc e Hot . This is to be contrasted with the GR-de 

Sitter solution: 3 Hq = U 0 => a(t) oc e,\I u °l it , which is 
obtained only for vacuum, p vac = Uq; p m = 0. 


The existence of the BD-de Sitter solution can be traced back 
to the work in the second item of Ref. [l8| . 
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The effective stiff-dilaton critical points (f2^ = 1): 

-P±:(a±,l) =>■ q± = 2 + V6a±, 



are also found, where the a± are defined in Eq. m- 


attractor solution, independent on the value of the BD 
coupling constant oj bd . 7 

In [17j | it has been shown that, in general, only for 
the specific exponential potential U{p) = M 2 exp <p, or 
for potentials which asymptote to M 2 exp <p, the ACDM 
model, which is associated with the GR-de Sitter point, 
can be an attractor of the Brans- Dicke theory. In 
other words, only BD cosmological models with U(ip) = 
M 2 expip, or U(<p) —> M 2 expip, can pass the present 
observational cosmological tests. 


C. The exponential potential 


IX. EMERGENCE OF THE GR-DE SITTER 
LOCAL ATTRACTOR 


In the general case when we have the exponential po¬ 
tential ®, since £ = 1 — A, is a constant, i. e., it 
can not be a phase space variable anymore, the plane- 
autonomous system of ODE-s simplifies to a single 
autonomous ODE: 


x = — 


(A + 2 + 2cj B d) x — a/§( 1 — A) 


1 + 2lu bd /3 

x(l + V6. 


x — LUbdX 


2 ) ■ (44) 

The critical points of the latter dynamical system are: 


\/273 (1 - A) 
A + 2 + 2cj bd 


X± Q!± 5 


(45) 


where the a± are given by Eq. (El). Notice that, since 
Xi ^ 0 (but for k = 1, in which case X\ = 0 and q = 
— 1), there are not critical points associated with constant 
ip = ipo. This means that the de Sitter phase with ip = 0 
(ip = const), U(ip) = const., i. e., the one which occurs in 
GR and which stands at the heart of the ACDM model, 
does not arise in the general case when A ^ 1. 

Hence, only in the particular case of the exponential 
potential (13U1) with A = 1 (£ = 0), which corresponds to 
the quadratic potential in terms of the JF-BD variables: 
V{<$) = M 2 cj) 2 , the GR-de Sitter phase is a critical point 
of the dynamical system (El) . In this case the critical 
points are (see Eq. (1751) 1: X\ = 0, x± = a±. Worth 
noticing that X\ = 0 corresponds to the GR-de Sitter 
solution 3 H 2 = M 2 exp ipo , meanwhile, the x± = a± , 
correspond to the stiff-fluid (kinetic energy) dominated 
phase: = 1. 

For small (linear) perturbations e = e(r) around the 
GR-de Sitter critical point: cc = 0 + e, e<C 1, one has 
that: e' = —3e => e(r) oc exp(—3r), so that it is an 


It is essential for the chameleon mechanism to work, 
that the chameleon potential V c h (0) given by Eq. (1171) . 
be a minimum at some Once this condition is met 
and the chameleon mechanism is in action, the BD field 
does oscillations around </>*, and it acquires an effective 
chameleonic mass as explained (see section ITlI AD . This is 
more easily seen in a cosmological context where, due to 
the expansion of the universe, the KG equation for the 
chameleonic BD field develops a friction term oc 4> which 
causes the oscillations to damp. 8 In consequence, the BD 
chameleon field does damped oscillations until it settles 
down in the minimum of the potential. Once this stable 
state is reached, the BD theory transmutes into general 
relativity with a cosmological A-term. 

The latter conclusion is true, however, only in the case 
when the density of the surrounding matter is a constant. 
Otherwise, if p = p(t) were a function of the cosmic time, 
we would not have a point of minimum of the chameleon 
potential, but a whole curve: ip* = ip*(p(t)). This would 
entail that 

^ °> 

the t-gradient of the dilaton will be aligned with the t— 
gradient of the density profile, so that we would not get 
general relativity but Brans- Dicke theory will be retained 
instead. 


' Around the stiff-matter solutions x = a± + e±: 

e±(r) oc e 3 ( 2 +'/®“±)' r j 

so that, if assume non-negative cj B d > 0, the points x± are always 
past attractors (unstable equilibrium points) since 2 + y/6a— > 
0. For negative cjbd < 0, these points are both past attractors 
whenever cj bd < — 3/2. In this latter case, for — 3/2 < cj bd < 0, 
the point x+ is a past attractor, while the point X— is a future 
attractor instead. 

8 One important example is the BD-KG equation in a cosmological 
Friedmann-Robertson-Walker space. 
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A. Transmutation of BD theory into GR 

In order to visualize how the transmutation of BD the¬ 
ory into GR-de Sitter theory operates, let us take a look 
at the Jordan frame BD action ©• We will assume that 
the given effective chameleon potential is a minimum 
at some constant expectation value of the JF-BD field 
cf> = (/>*, within a region of constant matter density. As 
commented above, these are necessary requirements for 
the transmutation of BD theory into general relativity. 

At the minimum <^>*, the gravitational sector of © 
transforms into the Einstein Hilbert action 


S* 


1 


cPxy/lgl (j)* 


R -2 


v(M 

(f>* 


(46) 


In a similar way, at cf> = </>*, the BD field equations 
©, ©, transform into the Einstein’s field equations of 
general relativity with a cosmological constant: 


G 


fllS 


JL T , ( m ) 

0. 


v {<!>*) 


9 


where, in the above equations, we have to set 


= Mi = 


8ttG 


a „d m =a, 


respectively. Hence, the chameleon effect in the BD the¬ 
ory warrants that the ACDM model is a local attractor of 
the Brans- Dicke cosmology. We have to point out, once 
again, that for this to happen it is indispensable that the 
given chameleon potential develops a minimum in a re¬ 
gion of constant density. Otherwise the former statement 
would not be true. This means the the GR-de Sitter so¬ 
lution can be, at most, a local attractor (it can be also 
a saddle point), since a global attractor would entail a 
constant density of the cosmological background. But, 
as known, any form of matter in a cosmological setting, 
but for the quantum vacuum with constant density p vac , 
would have, necessarily, an evolving density p = p(t) (see 
subsection IIXBI) . 

In order to compare with the results of the dynami¬ 
cal systems study (see section IVIIII) . it is easier to work 
with the dilatonic variable. Hence, assuming that the 
chameleon potential U c h(p) in Eq. (1271) . is a minimum at 
some constant p = p*, within a region of constant den¬ 
sity p , at the minimum the SF-BD equation of motion 
© decays into the Einstein’s GR equation: 

G^ = e~^T^-g^U(^), 

where e^* = M 2 L , and A = U(p*). In terms of the 
FRW metric we can say that the cosmological SF-BD 


equations (1M1) , transmute into the standard cosmological 
equations: 


3H = e~v 
w 


H = — 


Prr, 

+ 1 


+ U (<£>*) — Mp L p m + A, 


Pm = - - 


+ 1 


Mi 


Let us assume, for definiteness, that the spacetime 
dynamics is driven by the exponential potential 
U(p) = M 2 exp(Xp). Then, as shown, for A > 1, the 
corresponding chameleon potential m is a minimum at 
(see equation (l32]) J: 


p* = In 


P 

2M 2 (A — 1) 


1 

A+l 


Due to the expansion, i. e., to the friction term in 
the FRW-KG equation oc Hip, within the region with 
constant p the dilaton does damped oscillations around 
this p* until, eventually, the SF-BD theory gives way 
to Einstein’s GR with a cosmological constant. This re¬ 
sult is quite independent of the assumed dilaton potential 
U(p). In the discussed case (the exponential potential 
U(p) oc expyj), for instance, the result is independent of 
A (A > 1). 

How to reconcile this result with the result of section 
Em] (see Ref. 0,113)? In section IVIIII by means of the 
tools of the dynamical systems theory, it has been shown 
that, only for the exponential potential U(p) = M' 2 exp p 
(this corresponds to the quadratic potential V(<fi) oc <fr in 
terms of the JF-BD field), or for potentials that asymp¬ 
tote to M 2 exp p, the GR-de Sitter solution is an attrac¬ 
tor of the BD theory. Meanwhile, in this subsection it has 
been shown that the GR-de Sitter space is an attractor 
of SF-BD theory, provided that the effective chameleon 
potential is a minimum within a region of constant den¬ 
sity of matter. This result has been shown for the most 
general exponential potential U(p) oc exp(Ay>), indepen¬ 
dent of the free parameter A. I. e., the ACDM model 
is an atractor of SF-BD theory not only for the specific 
exponential potential U(p) oc expp, but for arbitrary 
exponentials. 

The answer to the above question is simple enough: in 
a cosmological setting, but for the quantum vacuum, the 
density of the background matter is an evolving function 
of the cosmic time t: p = p(t). Hence, assuming that the 
chameleon potential U c h(p) is a minimum at some 


P* = <P*(p(t)) => P* = P^ 0, 

dp 

which means, in turn, that at the minimum of the 
chameleon potential what we have is Brans-Dicke theory 
and not GR. Besides, since p * = p*{p(t)) is an implicit 
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function of the cosmic time, then the resulting effective 
mass squared m 2 ^ = m 2 (p(t)). would be so. Conse¬ 
quently, this can not be a stable state of BD theory, since 
in the latter the masses of particles (including any field 
excitations) are constants. This entails that the GR-de 
Sitter solution can be, at most, a saddle point, unless ei¬ 
ther the potential U(p) oc exp p, or it asymptotes to the 
exponential U{p) —> exp p, as it has been shown in sec¬ 
tion [VTTT] In general, only for regions of constant density 
of matter, for instance, regions filled with quantum vac¬ 
uum, the GR-de Sitter space can be a (local) attractor 
of the BD theory. 

B. Damped oscillations of the dilaton 

In order to illustrate our arguments above, let us inves¬ 
tigate the oscillations of the dilaton around the minimum 
of the exponential potential in a FRW spacetime with flat 
spatial sections (see section fVIIIl) . As before we shall as¬ 
sume that the FRW spacetime is filled with pressureless 
dust. The corresponding SF-KG equation © reads: 


is the effective mass squared of the perturbation, then, as 
long as p m is a function of the cosmic time, the effective 
mass is so: (t). But then, this mass would 

not be a standard mass of an elementary perturbation in 
SF-BD theory, since any physically meaningful mass pa¬ 
rameter in the SF/JF formulations of the Brans-Dicke 
theory should be independent of the spacetime point. 
Otherwise, as mentioned, coexistence of particles with 
constant mass and elementary perturbations with point- 
dependent mass, would entail serious problems with the 
equivalence principle. 

As a consequence of the above analysis, in a cosmo¬ 
logical setting, the chameleon effect makes sense only 
during small enough intervals of cosmic time, so that 
the background energy density may be considered as a 
constant. Only if the background were the quantum vac¬ 
uum with constant density p vac during the course of the 
cosmic evolution, the chameleon effect may have cosmo¬ 
logical sense. During a small time-interval St ~ 0(1), 
H(t) « H 0 ( 1 + a 0 5t), and p m (t) « p° m (l + faSt). Sub¬ 
stituting these equations back into the SF-KG equation 
above, and keeping the terms linear in the perturbations 
Sip ~ St ~ 0(1), one obtains: 


ip + 3 Hip + p 2 = 

2M 2 (A — 1) 
3 + 2w bd 



Pm e v 
2M 2 (X — 1) 


.(47) 


The second term in the LHS of this equation, i. e., the 
friction term oc Hp, is due to the implicit consideration of 
the curvature effects in FRW cosmological backgrounds. 

Next we expand this equation around the minimum of 
the chameleon potential (3D: 


Sip + 3 H 0 Sp « - (m° J“ Sp, 

2 (A + 1)[2M 2 (A-1)]^ t 

\ m vJ o I 9, . 


(49) 



Thanks to the friction term 3Hq Sip , the solution of the 
linear approximation to the SF-KG equation above, are 
the damped oscillations: 


e Xv ‘ - 


Pm & 


2M 2 (A — 1) 


= 0 


ip* = In 


_2M 2 (A — 1)_ 
Hence, in Eq. ED, we make the replacement: 


(48) 


p* + Sip —> p, Sp — > ip, Sp — > ip, 

where the perturbation Sp -C 1, i. e., Sp ~ Sp ~ Sp ~ 
0(1), etc. Up to linear terms, the corresponding KG 
equation for the perturbation reads: 


Sp + 3HSp = 


(A + 1) [2M 2 (A-1)] 
3 + 2w bd 


Pm 1 Sp. 


The problematic terms in the above equation are the 
time-dependent terms H(t ) and p m (t). Since the RHS 
term can be written as —m^*Sp, where 


(A + 1) [2M 2 (A-1)]~ 
3 + 2lo bd 


Sp(t) = Ce~~ 4 sin (w t + $), (50) 

where $ is an arbitrary phase which depends on the ini¬ 
tial conditions chosen, C is an integration constant, and 



This solution is valid only during small enough interval 
of the cosmic time St <C 2/3Hq ~ lifetime of the universe. 
Hence, the damped oscillations that transmute the SF- 
BD theory into Einstein’s GR with a cosmological con¬ 
stant, which are valid for any chameleon potential with a 
minimum, can not be an equilibrium configuration in the 
equivalent phase space. These can be, at most, a local 
attractor or, even, a saddle point. A global attractor in 
the phase space is a stable equilibrium point which does 
not depend on how far into the past the initial condi¬ 
tions for an specific orbit are given. As shown in section 
IVIII1 this is obtained only for the exponential potential 
U(p) oc exp p, or for potentials that asymptote to expy;. 
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Our conclusion is that one should trust the results 
of the study based on the dynamical systems theory 
(see section IVIIII and the references [IalI3)i where the 
chameleon effect is not explicitly considered, over the re¬ 
sults based on the above explained transmutation, since 
in a cosmological setting the latter effect is valid only 
during short stages of the cosmic evolution, where the 
background density of matter may be considered almost a 
constant value and, consequently, it can not give rise to a 
global attractor. This does not mean that the chameleon 
effect is not physically meaningful, but rather that, in 
a cosmological setting, its consequences are of limited 
reach. 

Our results here differ from those in Ref. 0 , where 
the cosmological chameleon is also investigated, but not 
within the context of the BD theory, and not in the 
Jordan/string frames, but in the Einstein frame. In 
that reference, the chameleon obviously violates the Ein¬ 
stein’s equivalence principle (EEP), since its effective 
mass evolves during the course of the cosmic evolution. 
This might be related with the fact that in 0 (see aslo 
1), the physically meaningful matter density is not the 
one measured by co-moving observers in the EF (this is 
not the conserved one in the EF), neither the conformal 
one, but a density which does not depend on the dilaton. 
In [ly|, the matter fields couple to the conformal metric 
g^l = exp(2/3i(f>/M PL ) g^, while the density of the non- 
relativistic fluid measured by EF co-moving observers is 
denoted by pi. It is assumed that what matters is the 
(/(-independent density pi = pi exp(3/3i</>/M PL ), which is 
the one conserved in the EF. While this choice may not 
be unique, in the Jordan as well as in the string frames of 
the Brans- Dicke theory one does not have this ambigu¬ 
ity: the matter density measured by JF(SF) co-moving 
observers p, is the one conserved in the Jordan/string 
frames and, additionally, does not depend on the BD- 
field (on the dilaton). 


X. THE COSMOLOGICAL CONSTANT 
PROBLEM 


The chameleon transmutation of BD theory into GR 
within regions of constant density of matter, poses a new 
problem that can be stated as follows. Let us imagine two 
separate spatial regions: one region with mean (constant) 
density of matter po, and another region with (also con¬ 
stant) density poo. Let the chameleon potential U c h(p) 
to have different minimums at each one of these regions: 
<po and <Poo, respectively. We have that 

1 1 


pVO — 

Po 

A+l 

pVoo _ 

Poo 


2M 2 (X — 1) 

J e 

2M 2 (X — 1) 


The resulting problem is that we will have two differ¬ 
ent values of the measured “Planck mass” M 2 L = exp ip. 
Their ratio can be written as: 


Ml 

(M°°)p L V/W 


(51) 


An also related problem is the well-known (old) cos¬ 
mological constant problem [2^, flbUlsl : why the present 
energy density of vacuum: p™™ ~ 10 -48 GeV 4 , is so 
small compared with the expected (theoretical) value: 
Pvac ~ 10 72 GeV 4 ? In our example the present (cosmo¬ 
logical) energy density will be identified with the mean 
density Poo, meanwhile, the theoretical expected value 
may be the mean density po- In consequence, we have 
that 


A o = U(<P o) = ( Po_\ X+1 ,„\ 

Aoo U(<p oo) VPoo/ 

where Ao = U(ipo) and A^ = U(tp oo), are the values of 
the cosmological constant measured in each one of the 
regions. 

Hence, we have two related problems that, for the ex¬ 
ponential potential case, are expressed by equations m 
and (15^1) . respectively, which are associated with the oc¬ 
currence of the chameleon effect in the SF-BD theory. A 
way out may be to assume large A > 1 (A > 119 to be 
precise). In this case we have 


( M °)l 

(M°°) 2 L 


= 1CV+ 1 


1, 


(53) 


i. e., the measured value of the Planck mass is almost 
the same in both regions. At the same time, there can be 
a huge difference between the value of the cosmological 
constant measured in one region when compared with the 
value measured in the other region: 


Ajl 120A i nn 

-2- = IOt+t ~ 10 120 , (54) 

Aoo 

as observed. In this case the origin of the cosmological 
constant problem may be attributed to the chameleon 
effect itself. 

Notice that, in order to get a consistent picture, here 
we have to set A 1, while the global GR-de Sitter 
attractor occurs for A = 1 (see section IVIIII) . 
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XI. JORDAN/STRING AND EINSTEIN 
FRAMES 

The conformal transformations controversy HH - fETi l is 
still a subject of debate [HH-HH. Under a conformal trans¬ 
formation of the metric: 


Let us assume, for simplicity, a single matter species, 
so that g$ = e p,p g^. The SF dilaton gravity action 
plus the matter fields reads: 


nip _ 


d i xsf\g\e' p R - w BD (dip) 2 - 2U(tp) 


9„u = =► v¥l = U" 4 v^j, (55) 


+2 J d A x\f\g\e 2 P v C m (x,dx,^ v g,iv) 


with SI 2 = </>, the JF DB theory, depicted by the action 
CD, is transformed into the Einstein frame formulation 
of BD theory: 


St = j d 4 xy/\f\ R~ (wb D + |j 

-2<T 2 U(0)] +2 J d^x^/Mt- 2 Cm, (56) 


where, as before, C rn = C m (x,c*x, ^ 1 9pv), is the La- 
grangian of the matter fields x, which are coupled to the 
JF metric g ^ = (j>~ l g 

There is an ongoing debate on whether theory © or 
theory HD, has the physical meaning [49l - l54j | . In this re¬ 
gard, the chameleon effect has been investigated, exclu¬ 
sively, in the Einstein frame version of the scalar-tensor 
theory. This fact leaves room for a certain suspicion on 
the real physical meaning of this effect. This is why we 
are discussing in the present paper the chameleon effect 
from the perspective of the Jordan and string frames for¬ 
mulation of Brans-Dicke theory. Although we shall not 
discuss here the very warped subject of the conformal 
transformation’s issue, nevertheless, certain features re¬ 
lated with these transformations will be discussed within 
the context of the chameleon effect. 


A. String—frame BD theory 

The string frame formulation of the dilaton-gravity ac¬ 
tion [55| admits, in principle, that different matter fields 
may couple to different conformal metrics, instead of the 
SF metric g^: 


d~-m ^X(i), SX(i) > g\xv^ 1 

where the matter fields X(i) couple to the conformal met¬ 
rics: 


« 

fits 


= e PiV g 


(57) 


The /3j-s are different coupling strengths of order unity. 


where we have dropped the subscript i, so that, as before, 
X stands collectively for the matter fields. For /3 = 0 we 
recover the string frame formulation of the BD theory 
given by the action ©. 

Under a conformal transformation of the metric (15511 
with n 2 = e^, the above action is transformed into the 
EF action: 


r-(u bd + -) (dv) 2 -2U(v) 


St = J d 4 x ^\= g I 

+2 f d 4 xy/\f\e 2( - p ~ 1)ip Cm (x)^X,e (/3_1)¥, fl/i^) >(59) 


where U(<p) = e ¥> [/(<^). Besides, under (l55l) (U 2 = e v ), 
the stress-energy tensor of the matter fields 


7 A m ) — 9_ 

± p,v ~ z 


VW\ eW-VvCm (X, Bx, e^-^gpv) 


\/M dg 


fils 


transforms like = e V T^\ where 


7 d m ) _ _o 


VWl{x,9x,g^ 


V\i\ d 9‘ 


fits 


The coupling of the matter fields to the conformal met¬ 
ric and not to the SF metric g means that the stress- 
energy is not separately conserved in the string frame, 
but that there is exchange of energy-momentum between 
the matter fields x and the dilaton: 


where T= g^T^ is the trace of the stress-energy 
tensor of matter. Under the conformal transformation 
HD with fl 2 = e the above “continuity” equation is 
transformed into: 


= P—ld v <fTW, 

so that the energy is not conserved in the Einstein frame 
neither. 
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Setting /3 = 0 in the former equations, amounts to 
the assumption that the matter fields are minimally cou¬ 
pled to the SF metric, which corresponds to the standard 
Brans-Dicke theory depicted by the dilatonic (SF) vari¬ 
ables. This means, in turn, that, in the string frame, the 
standard conservation equation is satisfied: = 0, 

while in the EF the stress-energy is not separately con¬ 
served: 


(155b . and in the related equations, one has to set 0 = 0 
as well. This means that we will be dealing with the 
EF-KG equation m- Following the same procedure of 
sections mm and IVIII we get that the effective mass 
squared of the dilaton in the Einstein frame is given by: 
mis = d 2 U ch . 

In the Einstein frame the chameleon potential is de¬ 
picted by: 


V LLV 



(60) 


E4h 


2 

3 + 2w bd 


U(<p) 


1 J'( m ) 

4 


(62) 


Instead = 0, where 


= ( w BD + g 


d„(pd„ip- i g„„ (dip) 2 


- 9n»U(ip), 


so that 


y Mj'(v’) = 


W BD + | ) VV - 9 V U 


d v ip. 


These equations lead to the EF-KG equation: 


V> = 


3 + 2o; b 


d v U + i 


(61) 


where there is subtlety that is to be explained. Assuming 
that in the Jordan/string frames the masses of particles 
mo, and the related quantities such as the density of mat¬ 
ter po, do not depend on spacetime position, their con¬ 
formal cousins in the Einstein frame fho(ip) = fI _1 mo = 
exp(— ip/2)mo, and po(ip) = O~ 4 po = exp(— 2<p)po, will 
be point-dependent quantities. Hence, for instance, 


d v Po = -2p 0 , d 2 p 0 = 4p 0 , ■ (63) 

etc. Then, since in the EF what the co-moving observers 
with 4-velocity = <5 q measure, is the energy density 

po = T^u" = T 0 ( 0 m) , 


which can be obtained, alternatively, by performing the 
conformal transformation (1551) , in Eq. 0, with Q 2 = e v . 

If choose 3 = 1 in the equations above, then the stress- 
energy of the matter fields is conserved in the Einstein 
frame, and the EF action (15U1) coincides with the stan¬ 
dard Einstein-Hilbert action of general relativity. In this 
case the action (l58l) will amount to the st ring frame (con¬ 
formal) formulation of general relativity [561 ]. 


and assuming, as before, that in the EF we have a pres¬ 
sureless (non-relativistic) fluid T= —po, the deriva¬ 
tive of the chameleon potential (l62l) : 


dm Uclci — 


3 -(- 2w b 


d<pU(<p) - \po(p) 


coincides with the RHS of Eq. flEED - 

Let us consider the SF exponential potential (13(J1) . so 
that, in the Einstein frame: 


XII. EFFECTIVE MASS OF THE 

CHAMELEONIC DILATON IN THE EINSTEIN 
FRAME 

Let us to compute the effective mass m^, of the 
chameleonic dilaton in the Einstein frame. This is the 
mass which determines the range of the Yukawa-like 
modification of the Newtonian potential 

e —r/A eff 

ALbv oc-, 

r 

where A e ff = is the Compton length of the dilaton 
in the Einstein frame. 

For definiteness we will consider the Brans-Dicke the¬ 
ory, which, in the string frame is given by ©, i. e, by 
the action (1581) with 0 = 0. Hence, in the EF action 


U(ip) = M 2 e {x ~ 1)tp . (64) 

The resulting EF chameleonic potential (16^1) . reads: 


t7ch 


2m2 r (a-i)c P , poM 
3 + 2w bd [ 4M 2 


(65) 


It is a minimum at the value ip * which solves the equation 

(A-l)cp. _ Po(P*) 

2M 2 (X — 1) ’ 

i. e., since po(p) = exp(— 2ip)p 0 , where po is a constant, 
then 


p* 


A ■ 


■ In 


Po 


2M 2 (X — 1) 
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This position of the minimum in the (^-direction coin¬ 
cides with the one in the string frame in Eq. (l32l) . We 
have that: 


Po(+*) = [4M 4 (A - l) 2 ] A+1 p 0 A+1 . (66) 

The effective chameleonic mass squared of the dilaton 
in the Einstein frame ?fi 2 ff = , is given by the follow¬ 

ing expression: 


= d l Uc h (^*) = 


(A + 1) p 0 {tp*) 
3 + 2w B d 


(67) 


where we prefer to keep the dependence on po, which is 
the density of matter measured by EF co—moving ob¬ 
servers, instead of on the SF constant pq. In other 
words, the effective EF mass squared of the dilaton de¬ 
pends linearly with the density of the surrounding matter 
ju 2 ^ oc poj as measured by EF co-moving observers. 

If desired, the EF mass squared ?u 2 ^ in Eq. m can 
be straightforwardly written in terms of the SF constant 
po as well, by substituting (l66l) in Eq. (l67l) : 


A. The SF and the EF are equivalent 
representations of the same theory 


According to one of the most popular points of view 
on the meaning of the conformal transformations of the 
metric [fol [5lj ]. both conformal frames: the SF/JF and 
the EF, are physically equivalent. The basic idea is that 
the two conformal frames are physically equivalent pro¬ 
vided that in the EF the units of time, length, mass, and 
derived quantities are allowed to scale with appropriate 
powers of the conformal factor f2. If one adheres to this 
point of view, as we temporarily do here, one must ac¬ 
cept that, instead of a system of units rigidly attached 
to the spacetime, the Einstein frame contains a system 
of units that depend on the spacetime point 5l| (for an 
alternative point of view on the equivalence see [531]). 

Let us consider that the mass of a given field, say, to/, 
is a constant in the string frame. Under a conformal 
transformation g =>- to/ = fU 4 TO/, so that, 
in the EF the mass of the field will be point-dependent, 
provided that in the SF the mass of the particle to/ is a 
constant. The fine point is that, in an experiment what 
one measures is the ratio fhf/fh u between the mass of 
the field and an arbitrarily chosen mass unit fh u . Hence, 
in the EF it is the conformal invariant ratio 


m l, = 


(A + 1) [4M 4 (A- l) 2 ]' 
3 + 2w B d 


■Po 


x-i 
+ i 


( 68 ) 


to / Q 1 rrif to / 
fh u f m u ’ 


(69) 


This equation confirms the transformation law for the 
mass under (1551) : fh v „ = fActually, under the 
conformal transformation (1551) with fl 2 = exp ip, 


— 2 —(p* 2 

= e m v>. 


Po 


2M 2 (A — 1) 


TO,, 


which, if substitute m 2 ^ from Eq. (1351) . yields (1551) . 


XIII. CONFORMAL FRAMES AND THE 
CHAMELEON EFFECT 


The controversy on which of the conformal frames 
where the Brans- Dicke theory can be formulated: the 
JF/SF or the EF, is physical, started with the paper 
by Dicke [49j, who stated that these are both equivalent 
representations of the same theory. However, no matter 
how trivial the question seems, the debate has not ceased 
for over the last 53 years [50l - l54l ] . In the present section 
we want to extend the debate to the discussion on the 
physical implications of the chameleon effect for the BD 
theory. 


what really matters & As a consequence, a measure¬ 
ment of the field’s mass yields a same value in the EF 
and in the SF. In other words, according to the present 
point of view, the given measurement can not differenti¬ 
ate between the different conformal frames. In order to 
illustrate this case, we take, for instance, the exponential 
potential U(p) oc exp(A(/>), and use the equation for the 
EF effective mass squared (1551) . instead of (1671) . In this 
case one has, 


/ A + 1 

|2M 2 (A — 1)J A+1 po (A+1J 

/ 3 + 2w B d 

H“ 1 ((/)*)to w 


/ A + 1 

V / 2M 2 (A — 1) 

*+T 2(A + 1) 

Po 

/ 3 + 2o; B d 

m u 



m v , 

m u 


In this case the chameleon effect is the same indepen¬ 
dent of the frame, and, provided that the observational 
bounds are met in one frame within a certain accuracy, 
these will be met in its conformal partner with the same 
accuracy as well. The contrary statement is also true. 
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B. The SF and the EF nest different theories 


C. The effective field-theoretical mass 


There is a radically different point of view on the 
(in) equivalence of the different conformal formulations 
of BD theory. According to this viewpoint, the SF-BD 
theory and the EF BD theory are totally different the¬ 
ories: i) same metric affinity (pseudo-Riemann spaces), 
ii) different equations of motion, i. e., different laws of 
gravity, iii) different experimental bounds to meet. In 
this case the conformal transformations of the metric are 
understood just as the deformation of one theory, that 
leads to the conformal one. 

In what regards to the chameleon effect within the BD 
theory, this viewpoint is supported by the equation (1671) 
for the EF dilaton’s effective mass squared. It is under¬ 
stood that, what really matters in the Einstein frame, 
are the quantities that EF co-moving observers with 4- 
velocity u 11 = Sq, measure. Co-moving observers in the 
EF see, for instance, exchange of energy-momentum be¬ 
tween the dilaton fluid and the standard matter fields, 
i. e., they “see” a fifth force which makes the motion of 
point-like particles to depart from geodesic motion in the 
EF metric g^. This is to be contrasted with what SF 
co-moving observers with 4-velocity = 5 q, see: the 
stress-energy tensor of each one of the matter species, is 
separately conserved. In other words, the SF co-moving 
observers do not see any fifth force effect. Besides, while 
the EF co-moving observers measure the energy density 
po = T 0 ( 0 m) = , the JF co-moving observers 

measure the constant energy density: po = Tp\ 

If assume that the present viewpoint is correct, i. e., 
that ffi v , oc y/po, then the EF-BD theory with expo¬ 
nential potential U(<p) oc exp[(A — l)<p], which is the EF 
cousin of U (<p) oc exp(Ap), improves the estimates of sec- 
tion lVIl bv several orders of magnitude. Actually, working 
as we did in section eh we obtain that 


m™ sm 


<p* 


I Pp 

^atm 

Po 


10 


-14 


(70) 


which is to be contrasted wit (E3l) . In the above equa¬ 
tion we have considered that p™* ~ 10 -31 g/cm 3 , and 
PQ tm ~ 1CV 3 g/cm 3 , are the critical energy density of the 
universe and the mean density of the Earth atmosphere, 
respectively, which are measured by EF co-moving ob¬ 
servers. 

In spite of the very discrete improvement of the EF- 
based estimate over the SF-based one for the exponential 
potential, it might represent a convenient observational 
signature, that may differentiate between the SF and EF 
formulations of the Brans- Dicke theory. 


According to Eq. (12511 . the effective field-theoretical 
mass squared of the dilaton in the string frame is given 

by 


'M = 


3 + 2uj b 


T ( m ) 

dp -d v U-—-e~‘<’ 


while, in the Einstein frame, given the chameleon poten¬ 
tial H, it is given by 


= dp ch = 


3 -\- 2uj b 


dp- T ("b 


where it has been taken into account that d v T t ' m> = 
-2 f( m ). There is no way in which one of the above 
equations can be obtained from the other one by means 
of a conformal transformation of the metric. Actually, if 
in the first equation above make the substitutions: 

u(ip) = n 2 u(p), r (m) = n 2 = e v , 


one obtains that 


m 2 p) 


2ev 

3 + 2uj bd 


dp + dp 


rp(m)- 

2 


(71) 


that can not be written in the form of the corresponding 
conformal transformation of a mass parameter: m 2 = 
D 2 m^, = e v m 2 . Only if the chameleon potential (in ei¬ 
ther frame) is a minimum at some <p = <p*, the obtained 
effective mass retains the sense of an ordinary mass pa¬ 
rameter. As a matter of fact, if U c h in (15^1) has a mini¬ 
mum, this means that, at the minimum 

dp ch = o => dp = 

so that, the equation ED can be written as: 


2e v 

m<p = 3 + 2w b 


dp - 


This means that the physical meaning is with the ef¬ 
fective mass which is obtained from second derivatives 
of a chameleon potential evaluated at the minimum. In 
contrast, the effective field-theoretical mass has limited 
physical implications. 


XIV. DISCUSSION AND CONCLUSION 

According to Eq. ED, for the exponential potential 
U(tp) oc exp(A<p), the effective mass squared of the dila¬ 
ton is given by: 
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m 


2 

<p* 


(A + 1) [2(A-1 )M 2 ] A + x 
3 + 2o; bd 


A 



(72) 


The dependence of the dilaton’s effective mass on the 
density of the environment ( p ), is the ba¬ 

sis for the chameleon effect. The interest in the latter 
as an effective screening mechanism, resides in that the 
dilaton field might have impact in the cosmological dy¬ 
namics and, yet, it might be effectively screened from 
the solar system and terrestrial experiments. This is in 
contrast to the mass due to an effective potential which 
does not depend on p. In this case, provided that the ef¬ 
fective mass is small enough, the dilaton may modify the 
cosmological dynamics. However, since the mass is small 
no matter where the experiments with the dilaton are 
performed, one has to care about the stringent bounds 
coming from the terrestrial fifth -force experiments, the 
solar system experiments to test the equivalence princi¬ 
ple, etc. 

In the case of the Brans-Dicke theory, one interesting 
application of the chameleon effect is related with the 
possible relaxation of the stringent lower bounds on the 
BD coupling parameter w BD > 40000. In order to illus¬ 
trate the way such a relaxation may occur, let us assume, 
as we did in section ED that the dilaton is immersed in 
the earth atmosphere with mean density p atm ~ ICG 3 
g/cm 3 . Assume, besides, that the 1 mm lower bound on 
the modification of the Newton’s law of gravitation ap¬ 
plies, i. e., that m~l ~ 1 mm, and, for sake of definite¬ 
ness, let us set A = 3. This value of the free parameter A, 
corresponds to the quartic potential V((f>) oc <p 4 in terms 
of the Jordan frame variables. If substitute the above 
numbers into Eq. (El 


m v» = 


4 (4M 2 ) 1/4 p 3 / 4 


3 + 2 uj-q 


m 1 [mm] 


0.1 \/3 + 2cJb 


2(4M 2 ) 1 / 8 (p[gr/ crn 3 ]^ 


3/8 ’ 


one gets that 


Wbd ~ lWM-1.5, (73) 

meaning that for M ~ 1, one can have that the BD 
theory with an w BD of order unity may be the correct 
description of gravity, without conflict with the existing 
experimental bounds, coming from terrestrial and solar 
system tests of the equivalence principle, fifth force and 
modifications of the Newtonian gravitational potential. 

The interesting question is whether the same theory 
can explain the observational data of cosmological origin 


as well. Worth noticing that the effective mass squared 
m 2 in Eq. CD). vanishes if A = 1, i. e., for the specific 
exponential U(ip) oc exp</>. This result is independent of 
the conformal frame: in the EF we have that (see Eq. 

USD), 


_ 2 (A + 1)[4M 4 (A-1) 2 ]^ 

m cp, = -X—W- P 7 ^, 

3 + 2w bd 

so that, the EF effective mass squared of the dilaton van¬ 
ishes if A = 1 as well. As a consequence, for A = 1, the 
chameleon effect does not arise in neither frame. This 
means that the dilaton mediated interaction of matter 
is of long-range (it is not screened), and we would have 
had detected the BD dilaton in terrestrial and solar sys¬ 
tem experiments. If, on the contrary, one assumes that 
A ^ 1, as we have shown above, the chameleon effect 
arises but, then the GR-de Sitter solution is not an equi¬ 
librium point of the dynamical system m, neither of 
the simpler ODE (l44l) . This means, in turn, that the 
ACDM model is not an attractor of the vacuum Brans- 
Dicke cosmology. 9 The absence of the GR-de Sitter crit¬ 
ical point entails that, in case it were an exact solution 
of the BD cosmological field equations (1511) . it were a 
non-generic unstable one, without interest for describing 
a long-lasting period of the cosmological evolution. 

Another essential aspect of the chameleon mecha¬ 
nism is related with the fact that a convenient effective 
chameleon potential U c h(ip, p), is a minimum at some </>*, 
where 


<9 v !7ch(v?*) = 0, m 2 ^ = <9 2 £/ ch ((/5*) > 0. 

In a cosmological setting this entails that the dilaton does 
damped oscillations around the expectation value ip *, un¬ 
til it stabilizes at the minimum of the potential. As a 
consequence, the Brans-Dicke theory with a chameleonic 
dilaton renders general relativity at the minimum of the 
potential. Notice, however, that, in general, the density 
of the environment p = p(t) is a function of the cosmic 
time. This is particularly true if one considers non-static 
distributions of matter. In this more general case what 
one has is not a point of minimum, but a whole curve 
ip*(t) = tp*(p(t)). Hence, the dilaton does not acquire an 
expectation value, but an expectation curve, whatever 
this means. In such a case, <p*(t) ^ 0 and one does not 
end up with general relativity but with Brans- Dicke the¬ 
ory instead. Hence, a necessary condition for the GR-de 
Sitter space to be a local attractor (or a saddle point) of 


9 As shown in the Ref. IH , this result is also true for Brans-Dicke 
cosmology with pressureless matter. 
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the BD cosmology is that the density of matter p be a 
constant within a given region. 

One gets the following big picture: Provided that 
the laws of cosmological expansion are governed by the 
Brans-Dicke theory with the potential U(ip), if there 
where regions IZt in our universe where the density of 
matter p,; were a constant (perhaps a different constant 
in each region), and assuming that within each one of 
such regions the potential U(ip) were a minimum at 
some <pi* = <pi*(pj) (different for each region also), then, 
thanks to the chameleon effect, after different periods 
of damped oscillations around the different minimums 
<Pi„, the BD theory transmutes into GR-de Sitter grav¬ 
ity within each one of these regions, with different sets 
of fundamental constants (M 2 Li ,Ai). Yet, if the poten¬ 
tial U (<p) were one that asymptotes to the exponential 
U(ip) — > exp ip, then the universe as a whole evolved into 
the final stable GR-de Sitter state. 

Now we are in position to put together the main re¬ 
sults of the present research, in order to get to physical 
conclusions (we shall include, also, results of |16l.ll7|). 

• It is a very difficult (perhaps hopeless) task to find a 
BD potential that allows to meet at once the terres¬ 
trial/solar system bounds and the bounds coming 
from cosmological considerations (section IVII) . 

• The potential U oc exp <p, is singular in the sense 
that the dilaton is strictly massless during the 
course of the entire cosmic history, since the corre¬ 
sponding Udi(tp) does not develop minimums (see 
Eq. (1721) with A = 1). Consequently, the chameleon 
effect does not arise. 

• Only for the specific exponential potential U(ip) oc 
exptp, or for potentials that asymptote to expy>, 
the GR-de Sitter solution is an attractor of the 
corresponding dynamical system EM3 (see the 
section IVIIII of the present paper). 

• In order to get a consistent local picture of trans¬ 
mutation of the Brans- Dicke theory into de Sitter- 
general relativity (see section IIXI) , assuming that 
the cosmological dynamics is dictated by the ex¬ 
ponential potential U(ip) oc exp(Ay>), it is required 
that X^> 1 (section DO). 
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XVI. APPENDIX: THE THIN-SHELL EFFECT 

The chameleon effect is primarily linked with the de¬ 
pendence of the effective mass of the BD scalar field upon 
the density of the surroundings. There is, however, an¬ 
other source of further screening of the fifth-force which 
can be associated with the BD scalar field. This is called 
as the thin-shell effect HI-0- This effect, which is sig¬ 
nificant only for large objects, is originated from the de¬ 
pendence of the strength of the coupling of the BD field 
to the matter, on the density of the environment as well. 
As a matter of fact, the BD scalar field outside such a 
large object, couples to matter with an effective coupling 
strength f3 e s which is much smaller than the “bare” cou¬ 
pling j3 (in the case of the BD theory considered in this 
paper /3 = 1). In order to consider the classical exam¬ 
ple [6], imagine an enough large spherical body of mass 
M c , homogeneous density p c and radius R c , which is im¬ 
mersed in an homogeneous environment of density poo- 
Denote by <p c and (f>oo, the values of the BD scalar field 
which minimize the chameleon potential I4h(/>) for p c 
and poo, respectively. The mass of the small perturba¬ 
tions of the BD scalar field about these minima will be 
denoted by m c and moo , respectively. We have that (see 
Eq. (El)): 

V c h(<j>(r)) = V eS (</>(r)) - p(r), 

3 + 2w bd 

where p = p c for r < R c , while p = poo for r > R c . 
The same for = m^(r) = yj 3|I4h(0(O) : m <t> = m c 
if r < R c , while m ^ = moo if r > R c . The resulting 
Helmholtz equation, 


The conclusion is that, if assume that the BD the¬ 
ory can explain the present acceleration of the expan¬ 
sion within the required accuracy, then the BD scalar 
field (our dilaton) has had to be detected in solar system 
and/or in terrestrial experiments. The contrary state¬ 
ment is also true: if the BD theory explains the local 
gravitational laws and, simultaneously, the dilaton is ef¬ 
fectively screened from terrestrial and solar system de¬ 
tection, then it can not explain the present accelerated 
stage of the cosmic evolution. 


1 d 2 d(<5</>) 
r 2 dr dr 


m| 6(j>, 


for the perturbations around <f> c and 4>oo, is subject to the 
following reasonable boundary conditions: i) the field is 
regular at the origin, and ii) far away from the body the 
field tends to its minimum value <f>oo, 


d(j>(r) 


= 0, lim <j)(r) = fioo 

r—> oo 


dr 


r =0 
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The solutions of the Helmholtz equation are 


where 


c/)(r < R c ) 


C\ sinh(?7i c r) 
r 


inside of the body, and 


C e~ rn °° r 

4>(r > R c ) = <f >oo H---, 

r 

outside it, respectively. Requiring, additionally, continu¬ 
ity of the field and its first derivative at the boundary 
r = R c , allows to fix the constants C\ and C 2 . In order 
to make estimates it is usually assumed that the density 
contrast between the body and the environment is large: 
4>c <j>oo, nic ^ rrioo. After the mentioned assumptions, 
the exterior solution above can be written as: 

M g-m^ir-Rc) 

4>(r > R c ) S3 2/3 e ff—--h 4>oo, 


3<p c = 9 

6 PcRl mlRV 

is the effective coupling of the chameleon field to the sur¬ 
rounding matter. For large objects with R c to” 1 , 
where to” 1 is the Compton wavelength of the chameleon 
field, the exterior solution is that of a point particle, but 
with a much weaker effective coupling to matter /3 e ff 1. 
This suppression mechanism is what is called as thin 
shell effect. The name is motivated by the fact that the 
BD (chameleon) field (j)(r ) is almost a constant through¬ 
out the bulk of the object except within a thin shell of 
thickness ~ to” 1 . Only that part of the body located 
within this thin shell, contributes to the fifth-force ex¬ 
erted on a test particle. This fact results in a very sup¬ 
pressed effective coupling /3 e g -C 1. 
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